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Quarkonium: an effective field theory approach 
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We discuss heavy quarkonium in the framework of pNRQCD, a QCD effective theory appro- 
priate for the study of heavy quark bound systems. 
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1 Introduction 

The experimental evidence that for heavy quark bound states like bb, cc, ... all the splittings are 
considerably less than the masses suggests that all the dynamical energy scales of these systems 
are small with respect to the quark masses. As a consequence the quark velocities v are small 
and these systems can be considered as non relativistic. The hierarchy of the scales is then 
the typical one of a non relativistic system. Called m the mass of the heavy quark, the quark 
momenta scale like mv and the quark energies like mv 2 . The inverse of the soft scale, mv, gives 
the size of the bound state, while the inverse of the ultrasoft scale, mv 2 , gives the typical time 
scale. However, in QCD another physically relevant scale has to be considered, namely the scale 
at which nonperturbative effects become important, which we will generically denote by Aqcd- 
Trickier is the situation in the gluon sector, but the binding interaction is essentially charac- 
terized by the same energy scale distribution. Therefore the dominant gluon interaction among 
heavy quarks appears "instantaneous". A potential picture should hold, at least in first ap- 
proximation, and the energy levels can be obtained by solving the corresponding Schrodinger 
equation. The study of the heavy quarkonium properties was initiated long ago, using potential 
models. A large variety of them exists in the literature and they have been on the whole phe- 
nomenologically quite successful. However, their connection with the QCD parameters is hidden, 
the scale at which they are defined is not clear, and they cannot be systematically improved. In 
spite of this, great progress has been achieved over the years by relating the potentials appearing 
in these models to some static Wilson loop operators^. This formulation is particularly suitable 



in QCD because it enables a direct lattice calculation and/or an analytic calculation in a QCD 
vacuum model S. However, there still remains the question of the extent this pure-potential 
picture is accurate. 

Here, we outline the kind of rigorous description of quarkonium that can be obtained in 
the framework of QCD effective field theories. The dynamical scales of the heavy quark bound 
system can be disentangled using effective theories of QCD. From the technical point of view, 
this simplifies considerably the calculation. From the conceptual point of view, this enables us to 
factorize the part of the interaction that we know from the low energy part which is dominated 
by nonperturbative physics. 

2 Potential NonRelativistic QCD (pNRQCD) 

The existence of a hierarchy of energy scales m S> mv 3> mv 2 , allows us to substitute scale by 
scale QCD with simpler but equivalent effective theories. Integrating out the hard scale m one 
passes from QCD to NonRelativisticQCD (NRQCD), while integrating out the soft scale mv 
one passes from NRQCD to potential NRQCD (pNRQCD) H. In this last effective theory only 
ultrasoft (US) degrees of freedom remain dynamical. The surviving fields are the Q-Q states 
(with US energy) and the US gluons. The Q-Q states can be decomposed into a singlet (S) and 
an octet (O) under color transformation. The relative coordinate r = xi — X2, whose typical 
size is the inverse of the soft scale, is explicit and can be considered as small with respect to the 
remaining (US) dynamical lengths in the system. Hence the gluon fields can be systematically 
expanded in r (multipole expansion). Therefore the pNRQCD Lagrangian is constructed not 
only order by order in 1/m, but also order by order in r. As a typical feature of an effective 
theory, all the non-analytic behavior in r is encoded in the matching coefficients, which can be 
interpreted as potential-like terms. 

The equivalence of pNRQCD to NRQCD, and hence to QCD, is enforced to the desired 
order in the multipole expansion by requiring the Green functions of both effective theories to 
be equal (matching). The matching between NRQCD and pNRQCD can be done order by order 
in I /ma. Moreover, if we restrict ourselves to the situation mv 3> Aqcd> we can, in addition, 
do the matching order by order in perturbation theory. 

The pNRQCD Lagrangian density is given at the next to leading order in the multipole 
expansion and at order (l/m)° by: 

£ P NRQCD = Trjst (id - V s (r) + . . .) S + 0* (iD - V a (r) + . . .) O j 

+gV A (.r)Tt {0 f r • E S + S f r • E o} + g^^-Tr {o f r • E O + f 0r • e} , (1) 

where R = (xi+X2)/2, S = S(r, R, t) and O = 0(r, R, t) are the singlet and octet wave functions 
respectively. All the gauge fields in Eq. (|l]) are evaluated in R and t. In particular E = E(R, t) 
and iDqO = idoO — <?L4o(R, t), O]. V s and V Q are the singlet and octet matching potentials 
respectively; Va and Vb are matching coefficients relevant at the next to leading order in the 
multipole expansion. 

We believe that pNRQCD is the most appropriate effective theory to describe the physics of 
heavy quark bound systems. Since pNRQCD has potential terms, it embraces potential models. 
Since in addition it still has ultrasoft gluons as dynamical degrees of freedom, it is able to 
describe non-potential effects. 

Indeed, the QED version of it, namely pNRQED, has been shown to correctly reproduce 
the non-potential effects that arise as 0(a 3 ) corrections to the binding energies of hydrogen-like 
atoms and positroniumi. Moreover, pNRQCD provides a new interpretation of the potentials 
that appear in the Schrodinger equation in terms of a modern effective field theory language. The 



potentials are nothing but r-dependent matching coefficients, which appear after integrating out 
scales of the order of the relative momentum in the bound state (~ mv ~ 1/r) or any remaining 
dynamical scale above mv 2 . Notice that, being matching coefficients, in the particular situation 
where no physical scales are between mv and mv 2 , the potentials generically depend on a 
subtraction point /j>. But, below the scale mv there are still dynamical US gluons, which are 
incorporated in the pNRQCD Lagrangian. The dynamics of gluons in pNRQCD is cut-off by 
/x, and this \x dependence will cancel the explicit /i dependence of the potential terms when 
calculating any physical quantity. 
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Table 1: Summary of the different kinematical situations. 

3 Results 

pNRQCD clarifies several long standing issues. In the case Aqcd < mv several situations have 
been examined! that we summarize in Tab. 1. Among them we quote: the proper treatment 
of the IR divergences in the singlet potential in perturbative QCD and the determination of the 
the leading log 3-loop contribution to iti; the construction of the suitable effective field theories 
for the ultrasoft degrees of freedom when mv S> Aqcd S> rnv 2 ; the model independent study of 
the short-distance limit of the hybrid static potentials. 

These results may become relevant in several situations. In particular, in accurate, model- 
independent, determinations of the bottom or top mass. From the perturbative point of view, 
the running of the singlet potential is the first step towards the full calculation of the leading 
log correction to the next-to-next-to-leading order results available at present for the above ob- 
servables. On the other hand our analysis provides a model-independent framework to estimate 
and parameterize nonperturbative effects. 

Finally let us briefly comment upon the situation where Aqcd ~ mv. Most of the ob- 
served charmonium and bottomonium states correspond to this situation, and hence the results 
presented here cannot be directly applied. This situation requires that the matching between 
NRQCD and pNRQCD is carried out nonperturbatively. Nevertheless, many of the features 
observed in this work survive in a nonperturbative analysis. The potential terms can still be 
obtained in an expansion in 1/m, while the multipole expansion can also be applied for the 
ultrasoft degrees of freedom. Moreover, in real QCD the Goldstone bosons (pions and kaons) 
remain dynamical at the ultrasoft scale producing non-potential effects. Work in this direction 
is in progress . 
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